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ABSTRACT 

The  response  of  stratospheric  climate  and  circulation  to  increasing  amounts  of  greenhouse  gases  (GHGs) 
and  ozone  recovery  in  the  twenty-first  century  is  analyzed  in  simulations  of  1 1  chemistry-climate  models  using 
near-identical  forcings  and  experimental  setup.  In  addition  to  an  overall  global  cooling  of  the  stratosphere  in 
the  simulations  (0.59  ±  0.07  K  decade-1  at  10  hPa),  ozone  recovery  causes  a  warming  of  the  Southern 
Hemisphere  polar  lower  stratosphere  in  summer  with  enhanced  cooling  above.  The  rate  of  warming  correlates 
with  the  rate  of  ozone  recovery  projected  by  the  models  and,  on  average,  changes  from  0.8  to  0.48  K  decade-1 
at  100  hPa  as  the  rate  of  recovery  declines  from  the  first  to  the  second  half  of  the  century.  In  the  winter 
northern  polar  lower  stratosphere  the  increased  radiative  cooling  from  the  growing  abundance  of  GHGs  is,  in 
most  models,  balanced  by  adiabatic  warming  from  stronger  polar  downwelling.  In  the  Antarctic  lower 
stratosphere  the  models  simulate  an  increase  in  low  temperature  extremes  required  for  polar  stratospheric 
cloud  (PSC)  formation,  but  the  positive  trend  is  decreasing  over  the  twenty-first  century  in  all  models.  In  the 
Arctic,  none  of  the  models  simulates  a  statistically  significant  increase  in  Arctic  PSCs  throughout  the  twenty- 
first  century.  The  subtropical  jets  accelerate  in  response  to  climate  change  and  the  ozone  recovery  produces 
a  westward  acceleration  of  the  lower-stratospheric  wind  over  the  Antarctic  during  summer,  though  this  response 
is  sensitive  to  the  rate  of  recovery  projected  by  the  models.  There  is  a  strengthening  of  the  Brewer-Dobson 
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circulation  throughout  the  depth  of  the  stratosphere,  which  reduces  the  mean  age  of  air  nearly  everywhere  at 
a  rate  of  about  0.05  yr  decade-1  in  those  models  with  this  diagnostic.  On  average,  the  annual  mean  tropical 
upwelling  in  the  lower  stratosphere  (—70  hPa)  increases  by  almost  2%  decade-1,  with  59%  of  this  trend 
forced  by  the  parameterized  orographic  gravity  wave  drag  in  the  models.  This  is  a  consequence  of  the  east¬ 
ward  acceleration  of  the  subtropical  jets,  which  increases  the  upward  flux  of  (parameterized)  momentum 
reaching  the  lower  stratosphere  in  these  latitudes. 


1.  Introduction 

This  paper  is  a  study  of  the  predicted  stratospheric 
climate  and  circulation  changes  in  a  multimodel  ensem¬ 
ble  of  chemistry-climate  model  (CCM)  simulations  of  the 
twenty-first  century.  The  simulations  were  carried  out 
in  support  of  the  2006  World  Meteorological  Organi¬ 
zation  (WMO)/United  Nations  Environment  Programme 
(UNEP)  Scientific  Assessment  of  Ozone  Depletion 
(WMO  2007).  Eyring  et  al.  (2006)  evaluated  the  ability  of 
the  CCMs  used  in  these  simulations  to  reproduce  the 
behavior  of  temperature,  trace  species,  and  ozone  in  the 
recent  past,  and  they  concluded  that  despite  a  huge 
range  of  skills  in  representing  the  different  processes 
assessed,  there  was  sufficient  agreement  between  the 
majority  of  the  CCMs  and  observations  that  some  confi¬ 
dence  could  be  placed  in  their  predictions.  Consequently, 
Eyring  et  al.  (2007)  were  able  to  analyze  projections  of 
the  ozone  recovery  in  the  simulations  of  the  twenty-first 
century,  but,  aside  from  noting  that  there  was  the  ex¬ 
pected  cooling  of  the  middle  atmosphere,  no  detailed 
analysis  was  carried  out  of  the  dynamical  changes  in  the 
simulations.  This  study  therefore  complements  that  of 
Eyring  et  al.  (2007). 

Climate  and  circulation  changes  in  the  stratosphere 
can  have  a  major  influence  on  the  ozone  layer  in  several 
ways.  Because  chemical  reactions  are  temperature  de¬ 
pendent,  a  systematic  cooling  of  the  upper  stratosphere 
can  slow  the  rate  of  photochemical  ozone  destruction  in 
that  region.  On  the  other  hand,  a  cooling  of  the  win¬ 
tertime  polar  lower  stratosphere  can  increase  the  occur¬ 
rence  of  polar  stratospheric  clouds  (PSCs)  and  enhance 
halogen-induced  ozone  loss  in  those  regions.  Circulation 
changes  can  have  an  impact  on  the  transport  of  ozone- 
depleting  substances  (ODSs)  into  the  stratosphere  and 
on  their  lifetimes  (Butchart  and  Scaife  2001),  in  addition 
to  affecting  temperatures  through  adiabatic  processes. 
Changes  in  both  transport  and  temperatures  will  af¬ 
fect  the  water  vapor  input  into  the  stratosphere,  while 
changes  in  dynamical  processes  will  affect  the  trans¬ 
port  of  ozone  within  the  stratosphere  and  across  the 
tropical  tropopause.  Therefore,  even  in  the  absence 
of  ODSs,  climate  change  can  alter  the  stratospheric 
ozone  distribution.  In  addition  there  is  now  evidence 
that  stratospheric  climate  and  circulation  changes  are 
a  factor  in  the  evolution  of  the  surface  climate  and  its 


variability  (e.g.,  Perlwitz  and  Harnik  2003;  Scaife  et  al. 
2005;  Baldwin  et  al.  2007). 

The  purpose  of  this  study  is  to  analyze  the  predicted 
stratospheric  climate  change  obtained  from  11  of  the 
CCMs  used  for  the  ozone  predictions  in  WMO  (2007). 
One  aim  is  to  establish  which  aspects  of  stratospheric 
climate  and  circulation  change  are  robust,  that  is,  model 
independent.  Although  all  of  the  models  are  CCMs,  the 
analysis  presented  in  this  paper  is  based  only  on  diag¬ 
nostics  that  are  obtained  from  the  variables  (e.g.,  winds 
and  temperature)  used  in  the  underlying  general  circu¬ 
lation  models  (GCMs),  with  the  exception  of  the  PSC 
area,  the  total  ozone,  and  an  age-of-air  tracer,  which, 
when  available,  is  used  to  diagnose  circulation  changes. 
In  contrast  to  previous  comparisons  of  stratospheric 
climate  predictions,  the  CCM  simulations  used  here  are 
all  transient  simulations  and  have  almost  identical  ex¬ 
perimental  setups  and  forcings,  though  there  are  dif¬ 
ferences  in  the  forcings  from  the  predicted  ozone  (e.g., 
Eyring  et  al.  2007)  and  other  radiatively  active  gases 
resulting  from  feedbacks  from  the  coupled  chemistry 
in  these  models.  The  analyses  update  or  extend  the  few 
other  published  intercomparisons  of  stratospheric  pre¬ 
dictions,  including  the  PSC  diagnostics  as  assessed  in 
Austin  et  al.  (2003),  the  temperature  changes  discussed 
in  Eyring  et  al.  (2007),  and  the  changes  in  tropical  up- 
welling  considered  by  Butchart  et  al.  (2006). 

A  brief  description  of  the  models  and  simulations  is 
presented  in  section  2.  Predicted  changes  in  the  strato¬ 
spheric  climate  (zonal  mean  temperature,  PSC  occur¬ 
rence,  and  eastward  wind)  can  be  found  in  section  3. 
Section  4  contains  an  analysis  of  the  circulation  changes 
with  the  downward  control  principle  of  Haynes  et  al. 
(1991)  used  to  diagnose  the  contributions  to  these  changes 
from  resolved  and  parameterized  subgrid-scale  waves. 
A  summary  and  concluding  remarks  appear  in  section  5. 

2.  Models  and  simulations 

Results  have  been  obtained  from  11  CCMs  that  per¬ 
formed  twenty-first  century  reference/sensitivity  simu¬ 
lations  for  WMO  (2007).  The  resolution  and  vertical 
domain  of  the  models  are  given  in  Table  1  along  with  the 
representation  of  subgrid-scale  momentum  deposition 
in  the  middle  atmosphere  from  nonorographic  gravity 
waves.  All  of  the  models  apart  from  the  University  of 
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Table  1.  Resolution,  vertical  domain,  and  nonorographic  gravity  wave  drag  parameterizations  used  in  the  models.  The  models  are 
listed  alphabetically  by  name.  For  gridpoint  models  the  horizontal  resolution  is  given  as  the  latitude  X  longitude  grid  spacing  (°).  For 
spectral  models  the  horizontal  resolution  is  given  as  the  triangular  truncation  of  the  spectral  domain,  with  the  equivalent  gridpoint 
resolution  in  brackets. 


Model 

Horizontal  resolution,  No.  of  levels,  top  level 

Nonorographic  GWD 

AMTRAC 

2°  X  2.5°, 

48  L,  0.0017  hPa 

Alexander  and  Dunkerton  (1999) 

CCSR/NIES 

T42  (2.8°  X  2.8°), 

34  L,  0.01  hPa 

Hines (1997) 

CMAM 

T32  (3.75°  X  3.75°), 

71  L,  0.0006  hPa 

Scinocca  (2003) 

E39C 

T30  (3.75°  X  3.75°), 

Enhanced  horizontal  diffusion 

39  L,  10  hPa 

in  the  two  uppermost  layers 

GEOS  CCM 

2°  X  2.5°, 

Based  on  Lindzen  (1981),  Holton  (1982), 

55  L,  0.01  hPa 

and  Garcia  and  Solomon  (1985) 

MAECHAM4CHEM 

T30  (3.75°  X  3.75°), 

39  L,  0.01  hPa 

Hines (1997) 

MRI 

T42  (2.8°  X  2.8°), 

68  L,  0.01  hPa 

Hines  (1997) 

SOCOL 

T30  (3.75°  X  3.75°), 

39  L,  0.01  hPa 

Hines (1997) 

ULAQ 

10°  X  22.5°, 

26  L,  0.04  hPa 

Rayleigh  friction  (Smith  and  Lyjak  1985) 

UMSLIMCAT 

2.5°  X  3.75°, 

Warner  and  McIntyre  (1996) 

64  L,  0.01  hPa 

and  Scaife  et  al.  (2002) 

WACCM  (version  3) 

4°  X  5°, 

Based  on  Lindzen  (1981),  Holton  (1982), 

66  L,  4.5  X  10~6  hPa 

and  Garcia  and  Solomon  (1985) 

L’Aquila  (ULAQ)  model  include  a  parameterization 
of  orographic  gravity  wave  drag  (GWD),  though  in  the 
Meteorological  Research  Institute  (MRI)  model  the  pa¬ 
rameterized  orographic  waves  do  not  propagate  into  the 
stratosphere  and  in  the  Atmospheric  Model  with  Trans¬ 
port  and  Chemistry  (AMTRAC)  the  drag  from  these 
waves  is  set  to  zero  above  42  hPa. 

Because  the  models  are  CCMs  they  include  both  gas- 
phase  chemistry  and  heterogeneous  chemistry  on  aero¬ 
sols  and  PSCs,  though  different  PSC,  dehydration,  and 
denitrification  schemes  are  used.  Except  for  the  Cana¬ 
dian  Middle  Atmosphere  Model  (CMAM),  which  did 
not  include  a  parameterization  for  nitric  acid  trihydrate 
(NAT)  PSCs  (see  Hitchcock  et  al.  2009,  p.  484b),  all  of 
the  models  calculate  solid  PSCs  based  on  simulated 
temperature,  nitric  acid,  and  water  vapor  concentrations 
using  the  Hanson  and  Mauersberger  (1988)  formulas. 
All  of  the  models  include  the  calculation  of  liquid  sulfate 
aerosols,  but  again  they  use  different  schemes.  In  addition, 
most  models  [i.e.,  Center  for  Climate  System  Research 
(CCSR)/National  Institute  of  Environmental  Studies 
(NIES),  CMAM,  Goddard  Earth  Observing  System 
(GEOS)  CCM,  Unified  Model  Single  Layer  Isentropic 
Model  of  Chemistry  And  Transport  (UMSLIMCAT), 
and  Whole-Atmosphere  Community  Climate  Model 
(WACCM)]  calculate  liquid  aerosols  from  a  supercooled 
ternary  solution  (STS)  using  approaches  similar  to,  for 


example,  Carslaw  et  al.  (1995)  or  Considine  et  al.  (2000). 
Further  details  of  the  models  can  be  found  in  Eyring  et  al. 
(2006,  their  Table  1  and  references  therein). 

The  reference  (REF2)  and  sensitivity  (SCN2)  simula¬ 
tions  broadly  follow  the  specifications  given  by  Eyring 
et  al.  (2005).  Briefly,  these  are 

•  surface  concentrations  of  greenhouse  gases  (GHGs) 
follow  the  Intergovernmental  Panel  on  Climate  Change 
(IPCC)  Special  Report  on  Emissions  (SRES)  GHG 
scenario  A1B  (medium;  see  Nakicenovic  and  Swart 
2000); 

•  sea  surface  temperatures  and  sea  ice  distributions 
(hereafter  referred  to  as  SSTs)  are,  if  available,  taken 
from  an  underlying  coupled  ocean-atmosphere  gen¬ 
eral  circulation  model  (OAGCM)  simulation  using 
the  same  GHG  scenario;  otherwise  the  SSTs  are  taken 
from  a  simulation  of  the  Met  Office’s  Hadley  Centre 
Global  Environmental  Model  version  1  (HadGEMl) 
OAGCM  (Johns  et  al.  2006),  again  using  the  A1B 
GHG  scenario;  and 

•  surface  halogens  are  prescribed  according  to  the  Ab 
scenario  of  WMO  (2003). 

The  REF2  simulations  include  anthropogenic  forcings  but 
not  natural  forcings,  such  as  variations  in  solar  irradiance 
or  big  volcanic  eruptions.  On  the  other  hand,  the  SCN2 
simulations  include  an  additional  forcing  representing 
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a  repetitive  solar  cycle  under  volcanically  clean  aerosol 
conditions.  Those  models  that  do  not  reproduce  a  quasi¬ 
biennial  oscillation  (QBO)  also  include  a  relaxation  forc¬ 
ing  toward  the  observed  QBO  (Giorgetta  and  Bengtsson 
1999)  for  the  SCN2  simulations. 

Table  2  summarizes  the  simulations  used  in  this  study 
and  any  differences  from  the  REF2/SCN2  specification. 
For  three  of  the  models  [E39C,  Middle  Atmosphere 
version  of  ECFIAM4  with  Chemistry  (MAECHAM4- 
CF1EM),  and  UMSLIMCAT],  the  simulation  only  ex¬ 
tended  over  the  first  two  decades  of  the  twenty-first 
century,  and  therefore  the  twenty-first  century  trends 
calculated  for  some  diagnostics  from  these  simulations 
are  often  not  statistically  significant  and  for  clarity  are 
excluded  from  some  of  the  figures.  Additionally,  results 
for  some  models  are  absent  from  a  particular  figure 
because  the  relevant  data  were  not  available  in  the  cen¬ 
tral  database  (e.g.,  gravity  wave  drag  diagnostics  were  not 
available  for  every  model;  see  section  4).  The  dynamical 
predictions  of  many  of  the  individual  simulations  have 
already  been  reported  in  the  literature.  References  to 
these  predictions  are  included  in  Table  2  and  at  appro¬ 
priate  points  throughout  this  paper.  Flere  we  concentrate 
on  using  the  multimodel  comparisons  to  investigate  the 
robustness  and/or  uncertainties  of  the  predictions. 

3.  Stratospheric  climate  change 

a.  Temperature  trends 

Because  the  surface  temperature  response  in  the 
simulations  used  in  this  study  is  strongly  constrained 
by  the  prescribed  SSTs,  the  main  focus  is  on  strato¬ 
spheric  change,  and  only  the  temperature  trends  above 
300  hPa  are  shown  in  this  section.  As  expected  from 
the  prescribed  increases  in  the  concentrations  of  the 
GFIGs,  all  of  the  models  predict  a  global  cooling  of  the 
annual  mean  stratosphere  for  the  twenty-first  century 
(Fig.  la).  Apart  from  the  results  for  the  E39C  model  and 
MAECHAM4CHEM  below  50  hPa,  these  trends  are 
statistically  significantly  different  from  zero  at  the  95% 
confidence  level,  where  a  lag-1  autocorrelation  is  used  to 
account  for  the  nonindependence  of  the  residual  values 
about  the  trend  line  (e.g.,  section  4  in  Santer  et  al.  2008). 
In  general  there  is  good  agreement  between  the  models 
in  the  magnitude  of  the  cooling,  particularly  in  the  lower 
and  middle  stratosphere,  and  all  of  them  predict  in¬ 
creased  cooling  with  height  in  the  stratosphere.  This  is 
consistent  with  the  trends  obtained  from  earlier  tran¬ 
sient  simulations  of  the  stratosphere-resolving  GCMs 
without  ozone  feedbacks  (e.g.,  Butchart  et  al.  2000), 
time-slice  experiments  with  ozone  feedbacks  (e.g., 
Jonsson  et  al.  2004),  and,  more  recently,  the  simulations 
performed  for  the  IPCC  Fourth  Assessment  Report 


(AR4;  e.g.,  Cordero  and  Forster  2006).  In  the  middle  to 
upper  stratosphere  (10-1  hPa)  the  models  predict  a  max¬ 
imum  annual  global  mean  cooling  between  0.63  (MRI) 
and  1.08  (ULAQ)  K  decade-1.  When  only  the  tropics 
(30°N-30°S)  are  considered,  the  transition  from  tropo¬ 
spheric  warming  to  stratospheric  cooling  becomes  much 
more  pronounced  and  occurs  at  a  slightly  higher  altitude 
(cf.  Figs.  la,b).  Again,  the  trends  are  nearly  all  statistically 
significantly  different  from  zero  at  the  95%  confidence 
level,  except  for  the  MAECFIAM4CF1EM  trends  below 
50  hPa.  In  the  tropics  the  trends  predicted  by  the  E39C 
model  are  statistically  significant  apart  from  between 
115  and  70  hPa.  For  the  five  models  (AMTRAC,  CCSR/ 
NIES,  CMAM,  GEOS  CCM,  and  MRI)  that  simulated 
the  whole  of  the  twenty-first  century,  there  was  very  little 
difference  between  global  and  equatorial  trends  calcu¬ 
lated  for  the  first  50  yr  from  those  for  the  full  100  yr 
shown  in  Fig.  1. 

In  general,  in  the  high  latitudes  (i.e.,  poleward  of  60°) 
there  were  large  seasonal  variations  in  the  temperature 
trends  predicted  by  the  models.  Figure  2  shows  the 
vertical  profile  of  the  trends  in  the  mean  temperatures 
poleward  of  60°  for  the  two  solstice  seasons  for  both 
hemispheres.  In  both  seasons  the  high-latitude  trends 
from  the  short  20-yr  simulations  of  the  E39C  model, 
MAECHAM4CHEM,  and  UMSLIMCAT  were  not  sta¬ 
tistically  significant,  except  in  the  Northern  Flemisphere 
upper  stratosphere  and  mesosphere  during  June-August 
for  the  MAECHAM4CHEM  and  UMSLIMCAT  (the 
E39C  model  does  not  include  the  upper  stratosphere  or 
mesosphere  in  its  domain).  Because  of  this  lack  of  sta¬ 
tistical  significance,  the  results  from  these  three  models 
were  not  included  in  the  multimodel  mean  trends  shown 
in  the  figure  although  the  individual  model  predictions 
are  shown. 

For  the  austral  summer  all  but  two  of  the  models  predict 
a  warming  of  the  polar  upper  troposphere/lower  strato¬ 
sphere  (Fig.  2a),  similar  to  that  found  by  Son  et  al.  (2008) 
and  is  almost  certainly  a  consequence  of  the  recovery  of 
ozone  over  Antarctica  in  the  simulations  (Eyring  et  al. 
2007).  Moreover,  the  spread  in  the  amplitude  of  the  pre¬ 
dicted  upper-tropospheric/lower-stratospheric  warming, 
from  about  0.3  K  decade-1  for  the  CMAM  ensemble  to 
up  to  1  K  decade- 1  for  the  WACCM  ensemble,  appears  to 
be  directly  related  to  the  predicted  rate  of  ozone  recovery. 
Figure  3  shows  that  poleward  of  60°S  during  December- 
February  there  is  a  good  correlation  across  the  multi¬ 
model  ensemble  between  the  predicted  temperature 
trends  at  100  hPa  and  the  projected  trends  in  total 
ozone.  For  the  eight  models  that  simulated  the  period 
2000-50  the  correlation  coefficient  between  the  ozone 
and  temperature  trends  was  0.83  with  the  slope  of  the 
least  squares  linear  fit  (solid  line  in  Fig.  3),  implying 
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Fig.  1.  Vertical  profiles  of  the  annual  mean  (a)  global  and  (b)  tropical  (30°N-30°S)  temperature  trends  predicted 
for  the  twenty-first  century.  The  trends  were  calculated  from  a  least  squares  linear  fit  to  the  modeled  temperature 
time  series  for  the  years  of  2000-99  (solid  curves)  or  to  the  end  of  the  simulation  if  that  was  before  2099;  that  is,  2050 
for  the  SOCOL,  ULAQ,  and  WACCM  simulations  (dashed  curves)  and  2019  for  the  E39C,  MAECHAM4CHEM, 
and  UMSLIMCAT  (dotted  curves)  simulations;  see  Table  2.  The  extension  “_e”  to  a  model  short  name  indicates  that 
the  ensemble  mean  of  more  than  one  simulation  was  used  to  calculate  the  trend  for  that  model.  The  black  line  is  the 
multimodel  mean  with  the  gray  shading  denoting  ±  two  standard  error  (i.e.,  the  multimodel  standard  deviation 
divided  by  s/n ,  where  n  is  the  number  of  models). 

a  0.07  K  decade-1  increase  in  the  rate  of  warming  for 
each  1  DU  decade-1  increase  in  the  rate  of  ozone  re¬ 
covery.  Extrapolating  the  linear  fit  to  zero  ozone  trend 
suggests  a  warming  trend  at  this  level  of  0.13  K  decade-1 
attributable  to  the  GHG-induced  climate  change  (i.e., 

0.13  K  decade-1  is  the  intercept  in  Fig.  3  for  a  zero- 
column  ozone  trend).  Including  results  from  the  2000-19 
simulations  of  the  E39C  model  and  UMSLIMCAT  re¬ 
duces  the  correlation  between  the  ozone  and  tempera¬ 
ture  trends  only  slightly  from  0.83  to  0.81  and  has  very 
little  impact  on  the  slope  or  position  of  the  least  squares 
linear  fit  line.  Five  models  (AMTRAC,  CCSR/NIES, 

CMAM,  GEOS  CCM,  and  MRI)  also  provided  predicted 
trends  for  2050-99.  Again,  Fig.  3  shows  that  there  is  a 
good  correlation  between  the  models  (0.99;  but  note  the 
small  sample  size)  with  the  linear  fit  being  very  similar  to 
that  obtained  for  the  previous  50  years  (a  slope  of  0.06 
and  an  intercept  of  0.17  for  2050-99  compared  to  0.07 
and  0.13,  respectively,  for  2000-50).  Also,  Fig.  3  shows 
that  during  the  second  half  of  the  century  four  out  of  the 
five  models  project  a  slower  rate  of  ozone  recovery  (or 
increase),  and  hence  warming  at  100  hPa,  compared  to 
the  first  half  of  the  century.  On  average,  the  rate  of  ozone 


recovery  is  about  30%  less,  with  a  corresponding  decrease 
in  the  rate  of  warming  from  0.8  to  0.48  K  decade-1. 

The  lower-stratospheric  warming  is  accompanied  by 
cooling  of  the  middle  stratosphere  of  the  order  of  0.8 
K  decade  - 1  with  generally  good  agreement  between  the 
models  (Fig.  2a).  The  0.8  K  decade-1  cooling  of  the 
middle  stratosphere  is  more  than  would  be  expected 
from  the  changes  in  the  GFIG  concentrations  [cf.,  Fig.  la 
for  the  global  annual  mean;  see  also  Fig.  4b  of  Hitchcock 
et  al.  (2009)],  but  the  additional  cooling  accompanying 
the  warming  below  is  consistent  with  the  temperature 
response  that  Kiehl  et  al.  (1988)  and  Mahlman  et  al. 
(1994)  found  for  a  specified  ozone  hole.  An  explanation 
for  this  response  was  subsequently  provided  by  Manzini 
et  al.  (2003)  in  terms  of  enhanced  planetary  wave  drag  in 
December-February  from  the  delayed  breakdown  of 
the  vortex.  For  the  austral  winter  all  of  the  models  apart 
from  the  UMSLIMCAT  and  one  E39C  simulation  in¬ 
dicate  a  cooling  of  the  stratosphere,  which  increases  with 
height,  as  does  the  spread  between  models  (Fig.  2b). 

For  the  boreal  winter  there  is  good  agreement  be¬ 
tween  six  of  the  eight  simulations  that  extended  to 
2050  and  beyond,  with  almost  zero  trend  in  the  lower 
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Fig.  2.  As  in  Fig.  1,  but  for  the  temperatures  averaged  between  60°  and  90°S  for  (a)  December-February,  and 
(b)  June-August,  and  averaged  between  60°  and  90°N  for  (c)  December-February,  and  (d)  June-August. 


stratosphere  and  cooling  in  the  middle  and  upper  strato¬ 
sphere  (Fig.  2c).  The  MRI  and  ULAQ  models  are  the 
exceptions.  In  the  other  six  models  the  near-zero  future 
trend  in  the  lower  stratosphere,  which  also  extends  into 
spring  (not  shown),  contrasts  with  the  results  of  Austin 
et  al.  (2009)  who  found  that  when  the  same  models  were 


used  to  stimulate  the  recent  past  there  was  cooling  in  this 
region  in  February-April.  Although  this  difference  in  the 
past  and  future  responses  in  spring  suggests  that  the  pro¬ 
jected  recovery  of  Arctic  ozone  in  the  twenty-first  century 
(Eyring  et  al.  2007)  is  perhaps  having  a  direct  effect  on 
the  model  temperatures,  this  would  be  unlikely  during 
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Column  ozone  trend  [Du/decade] 


Fig.  3.  Scatter  diagram  of  the  trend  in  the  December-February  mean  column  ozone  at  60°- 
90°S  against  the  trend  in  the  100-hPa  temperatures  averaged  over  the  same  3  months  and 
region.  The  trends  for  2000-50  (crosses;  2000-19  for  E39C  and  UMSLIMCAT)  and  the  trends 
for  2050-99  for  the  five  models  for  which  the  simulations  extended  to  2099  (diamond  symbols 
with  model  short  names  in  italics).  The  solid  (dotted)  line  is  the  least  squares  linear  fit  to  the 
results  for  2000-50  (2050-99).  When  results  for  more  than  one  simulation  were  available  for 
a  particular  model,  the  trends  were  calculated  from  the  ensemble  mean. 


the  polar  night  (i.e.,  midwinter).  A  more  likely  expla¬ 
nation  for  the  winter  period  is  that  the  climate  change 
in  the  simulations  increases  adiabatic  warming  of  the 
Northern  Hemisphere  polar  stratosphere  so  that  it  more 
or  less  balances  the  extra  radiative  cooling  from  the 
increased  amounts  of  GHGs.  Consistent  with  this,  the 
mean  downwelling  in  the  lower  stratosphere  (—50  hPa) 
poleward  of  60°N  increased  by  0.6%  decade-1  for  the 
AMTRAC  and  CCSR/NIES  model,  0.8%  decade-1  for 
the  WACCM,  1.6%  decade-1  for  the  CMAM,  1.9% 
decade-1  for  the  GEOS  CCM,  and  3.5%  decade-1  for 
the  Solar  Climate  Ozone  Links  (SOCOL)  model,  though 
apart  from  the  CMAM  and  GEOS  CCM  results  these 
trends  are  not  statistically  significant.  In  the  other  two 
long  simulations  (MRI  and  ULAQ)  the  polar  downwell¬ 
ing  decreased  (and  also  was  not  statistically  significant), 
again  consistent  with  the  temperature  trends  in  Fig.  2c  for 
those  two  models.  However,  because  the  ULAQ  simu¬ 
lation  lacks  any  long-term  trend  in  its  surface  forcing  (see 
Table  2),  the  ULAQ  temperature  and  downwelling  trends 
are  less  certain  than  those  obtained  from  the  other  models. 
Therefore,  based  on  the  multimodel  ensemble,  the  in¬ 
creased  adiabatic  warming  of  the  Northern  Hemisphere 
polar  lower  winter  stratosphere  can  almost  certainly  be 
considered  as  a  robust  feature  of  climate  change.  For  the 
boreal  summer  there  is  reasonably  good  agreement  be¬ 
tween  the  models  (Fig.  2d)  and  the  trends  have  broadly  the 
same  features  as  those  for  the  global  mean  temperatures. 


b.  Polar  stratospheric  cloud  formation 

During  times  of  elevated  amounts  of  halogens,  changes 
in  stratospheric  temperatures  are  expected  to  have  a 
large  impact  on  polar  ozone  loss  through  the  formation 
and  occurrence  of  PSCs.  Following  Pawson  et  al.  (1999) 
and  Austin  et  al.  (2003),  Fig.  4  shows  a  time  series  of 
seasonally  accumulated  areas  for  the  Antarctic  (July- 
December)  and  the  Arctic  (November-April),  where  the 
daily  temperature  poleward  of  60°  is  below  195  K  (ap¬ 
proximate  NAT  and  STS  threshold  temperature,  upper 
row)  and  below  188  K  (approximate  ice  threshold  tem¬ 
perature,  middle  row)  for  the  models  (colored  lines)  and 
40-yr  European  Centre  for  Medium-Range  Weather 
Forecasts  (ECMWF)  Re-Analysis  (ERA-40;  black  dashed 
line).  The  quantity  used  in  Fig.  4  is  the  percentage  of  the 
horizontal  area  of  the  hemisphere  where  the  50-hPa 
temperatures  poleward  of  60°  are  below  the  NAT  and 
ice  PSC  formation  thresholds.  These  daily  percentage 
areas  are  then  accumulated  over  the  course  of  the  winter 
and  spring  (July-December  in  the  SH,  and  December- 
April  in  the  NH)  to  provide,  for  that  year,  an  estimate  of 
the  total  amount  of  NAT  and  ice  PSCs  in  units  of  per¬ 
cent  days,  as  shown  in  the  figure. 

In  the  Antarctic,  the  magnitude  of  the  multimodel 
mean  temperature-derived  NAT  (/’mat)  and  ice  (7’ice) 
PSC  area  agrees  well  with  the  ERA-40  data,  although 
it  should  be  cautioned  that  these  quantities  and  their 
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Fig.  4.  Seasonal  accumulation  of  the  (left)  Antarctic  (July-December)  and  (right)  Arctic  (November-April)  area 
at  50  hPa,  where  daily  temperatures  are  (top)  below  195  K  (approximate  threshold  temperature  for  NAT  and  STS 
formation)  and  (middle)  below  188  K  (approximate  threshold  temperature  for  ice  formation)  for  the  models  (colors) 
and  ERA-40  (black-dashed  line),  (bottom)  The  actual  PSC  area  in  the  models,  which  was  available  only  from  a  subset 
of  models  (CCSR/NIES,  MRI.  ULAQ,  UMSLIMCAT,  and  WACCM).  The  solid  black  line  is  the  multimodel  mean 
with  the  gray  shading  denoting  two  standard  error  as  in  Fig.  1.  In  the  lower  panel  UMSLIMCAT  and  WACCM  are 
excluded  from  the  multimodel  mean,  because  they  only  cover  a  short  period  and  produce  jumps  in  the  multimodel 
mean  if  included. 


trends  are  hard  to  derive  from  global  assimilation  data 
and  often  depend  on  which  analysis  or  reanalysis  data¬ 
sets  are  used  (e.g.,  Manney  et  al.  2003,  2005).  The  ma¬ 
jority  of  the  models  show  an  increase  in  7NAT  in  all  time 


periods  (the  entire  REF2  simulation,  1980-99,  2000-49, 
and  2050-99;  see  also  Table  3).  Only  the  UMSLIMCAT 
and  the  CCSR/NIES  model  simulate  negative  trends, 
but  these  negative  trends  are  not  statistically  significant 
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Table  3.  Linear  trends  in  Arctic  and  Antarctic  PSC  areas  at  50  hPa  for  each  model  derived  for  different  periods  (entire  time  period  of 
REF2  simulations,  1980-99, 2000-49,  and  2050-99).  Trends  that  are  statistically  significant  at  the  95%  confidence  level  are  shown  in  bold. 
In  addition  to  the  individual  models,  the  multimodel  mean  is  given. 


Model 

Period 

NH  60°-90°N  at  50  hPa 

SH  60°-90°S  at  50  hPa 

7n  AT /A  PSC 

7)ce  Actual  PSC  area 

7nat/APsc 

Tice  Actual  PSC  area 

ERA-40 

1979-2006 

-1.78/0.16 

0.68 

— 

-0.36/0.08 

4.87 

— 

AMTRAC 

1990-2099 

0.32/0.02 

0.0 

— 

0.72/0.02 

0.79 

— 

1990-99 

0.76/0.01 

0.0 

— 

2.18/— 0.04 

0.48 

— 

2000-49 

0.89/0.03 

0.0 

— 

0.66/0.02 

0.83 

— 

2050-99 

0.05/-0.02 

0.01 

— 

0.51/0.01 

0.58 

— 

CCSR/NIES 

1980-2050 

0.68/0.02 

0.03 

1.03 

-0.01/0.0 

0.35 

-0.29 

1980-99 

5.95/0.13 

0.03 

3.96 

1.01/— 0.02 

2.80 

1.85 

2000-49 

0.49/0.01 

0.01 

1.01 

0.35/— 0.01 

1.19 

0.00 

2050-99 

0.98/0.02 

0.01 

1.26 

0.26/0.00 

-0.06 

-0.01 

CMAM 

1960-2099 

-0.03/0.00 

0.0 

— 

1.33/0.04 

1.49 

— 

1980-99 

-0.54/ -0.06 

-0.08 

— 

4.01/0.12 

4.02 

— 

2000-49 

-0.53/-0.01 

-0.03 

— 

1.27/0.05 

1.46 

— 

2050-99 

-0.60/ -0.03 

— 

— 

0.11/0.00 

0.11 

— 

MRI 

1980-2099 

2.19/0.05 

0.79 

2.48 

0.67/0.01 

0.64 

0.75 

1980-99 

1.70/— 0.11 

-0.45 

1.85 

6.16/0.15 

4.4 

6.50 

2000-49 

2.42/0.05 

1.06 

2.62 

0.68/0.01 

0.73 

0.79 

2050-99 

0.29/0.09 

-0.02 

1.33 

0.78/0.00 

0.61 

1.10 

SOCOL 

1980-2050 

0.64/0.06 

-0.02 

— 

1.46/0.04 

1.77 

— 

1980-99 

1.22/0.16 

-0.15 

— 

6.41/0.15 

5.68 

— 

2000-49 

0.79/0.09 

-0.02 

— 

1.05/— 0.003 

1.47 

— 

ULAQ 

1980-2050 

0.63/0.03 

0.09 

0.70 

3.59/0.12 

3.48 

2.46 

1980-99 

0.18/0.01 

0.01 

0.14 

5.24/0.11 

7.51 

3.06 

2000-49 

1.58/0.07 

0.22 

1.77 

0.64/0.00 

0.52 

-0.62 

UMSLIMCAT 

1980-2020 

-0.35/0.04 

-0.15 

-1.41 

0.58/0.01 

0.80 

-1.10 

1980-99 

-1.3/— 0.14 

-0.47 

-2.4 

-4.03/0.04 

-3.59 

-3.33 

WACCM 

1980-2050 

-0.09/0.01 

— 

-0.01 

0.69/0.02 

— 

0.44 

2009-50 

-0.09/— 0.01 

— 

-0.01 

0.07/0.02 

— 

0.44 

Multimodel  mean 

1980-2099 

0.51/0.02 

0.0 

2.70 

1.43/0.04 

1.89 

0.93 

1980-99 

0.32/0.00 

-0.23 

0.75 

0.90/— 0.04 

2.01 

3.47 

2000-49 

0.08/0.00 

0.23 

0.65 

1.90/0.03 

1.49 

-1.35 

2050-99 

0.13/0.01 

0.67 

2.76 

0.19/0.01 

0.63 

1.20 

in  any  of  the  time  periods.  For  each  model  that  shows 
positive  trends,  the  increase  in  7NAT  is  largest  in  the 
1980-99  period.  A  further,  but  smaller,  increase  is  sim¬ 
ulated  in  the  first  half  of  the  twenty-first  century  and  an 
even  smaller  positive  trend  occurs  in  the  second  half  of 
the  twenty-first  century.  The  change  in  magnitude  of  the 
positive  trend  for  the  different  time  periods  is  consistent 
with  the  radiative  effects  of  simulated  ozone  depletion 
between  1980  and  1999  and  ozone  recovery  in  the  twenty- 
first  century,  in  addition  to  continued  cooling  of  the 
stratosphere  resulting  from  the  prescribed  increases  in 
GFIG  amounts  (Eyring  et  al.  2006, 2007).  A  similar  result 
is  found  if,  instead  of  the  percentage  of  the  hemisphere 
covered  times  their  duration  (days),  the  temperature- 
derived  NAT  PSC  area  is  plotted  (km2;  APsc)  at  50  hPa 
and  the  long-term  trend  is  calculated  from  the  maximum 
values  of  A  PSC  for  5-yr  intervals,  as  in  Rex  et  al.  (2006; 
see  also  our  Fig.  5).  Flitchcock  et  al.  (2009)  analyzed 
Apse  at  50  hPa  in  the  three  CMAM  simulations  used 
here  and  found  that  for  that  model  the  trends  in  APSC 


were  a  direct  radiative  consequence  of  changes  in  ozone 
and  GFIG  amounts  and  were  not  the  result  of  adiabatic 
changes  in  temperature.  The  trends  in  Tice  in  the  models 
show  similar  features.  Overall,  there  is  clearly  an  increase 
in  colder  winters  in  the  Antarctic  with  temperature  ex¬ 
tremes  below  195  and  188  K  favoring  NAT  and  ice  for¬ 
mation  in  the  years  of  ozone  depletion  (1980-99); 
however,  the  positive  trend  is  decreasing  over  the 
twenty-first  century  in  all  of  the  models  (see  Table  3). 

In  the  Arctic,  the  situation  is  not  as  obvious  as  in  the 
Antarctic.  There  are  large  differences  among  the  simu¬ 
lations  in  both  the  absolute  values  and  variability.  For 
example,  the  MRI  model  overestimates  ?nat  and  Tice 
compared  to  the  ERA-40  data.  This  result  is  consistent 
with  a  cold  bias  in  this  model  in  the  Arctic  winter  at  this 
altitude  (e.g.,  Fig.  1  in  Eyring  et  al.  2006).  Consistently, 
the  CCSR/NIES  model  is  also  on  the  high  side.  On 
average,  most  other  models  underestimate  7’N/\T  com¬ 
pared  to  ERA-40  and  the  multimodel  mean  for  7’NAT  is 
slightly  lower  than  that  derived  from  the  reanalysis  data. 
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The  models  are  more  consistent  in  the  simulated  trends. 
For  the  periods  of  1980-99,  2000-49,  and  2050-99,  the 
simulated  trends  in  Tnat  are  only  statistically  significant 
in  two  of  the  models  (AMTRAC  and  ULAQ).  Three 
models  (CMAM,  UMSLIMCAT,  and  WACCM)  simu¬ 
late  negative  trends,  but  again  none  of  these  trends  is 
statistically  significant  in  any  time  period.  In  the  CMAM 
the  lower  trends  in  the  last  decades  of  the  twentieth 
century  could  be  due  to  weaker  ozone  depletion  simu¬ 
lated  in  that  model  compared  to  the  observations 
(Hitchcock  et  al.  2009).  Consistently,  the  UMSLIMCAT 
and  WACCM,  which  also  simulate  smaller  ozone  de¬ 
pletion  than  that  observed  (Eyring  et  al.  2006,  their 
Fig.  15),  show  negative  trends  in  rNAT.  A  possible  rea¬ 
son  for  the  underestimation  of  ozone  depletion  by  the 
CMAM  is  the  absence  of  a  NAT  PSC  parameterization 
(Hitchcock  et  al.  2009).  All  of  the  other  models  simulate 
positive  trends  that  are  statistically  significant  only  in 
some  simulations  and  for  some  time  periods  (2000-49 
in  the  AMTRAC  and  ULAQ  model,  1990-2099  in  the 
AMTRAC,  and  1980-2050  in  the  CCSR/NIES,  SOCOL 
and  ULAQ  models).  The  multimodel  mean  indicates 
a  small  but  statistically  significant  increase  in  7  NAT  for 
the  second  half  of  the  twenty-first  century  and  a  statis¬ 
tically  significant  increase  of  0.51%  days  in  the  season¬ 
ally  derived  PSC  area  between  1980  and  2099.  The 
multimodel  mean  for  Tice  shows  an  increase  in  trends 
from  1980-99  to  2000-49,  which  is  mainly  due  to  the 
increase  in  7’ice  in  the  MRI  model.  However,  the  mul¬ 
timodel  mean  trends  that  are  given  in  Table  3  should  be 
carefully  interpreted,  in  particular  for  the  periods  of 
2050-99  and  1980-2099,  because  the  simulations  from 
different  models  cover  different  periods,  which  pro¬ 
duces  jumps  in  the  multimodel  mean,  for  example,  in 
2050.  The  trends  calculated  from  the  maximum  values  of 
APSC  for  5-yr  intervals  in  the  models  (Fig.  5,  lower 
panel)  are  all  smaller  than  the  trends  in  7NAT,  but  they 
agree  well  in  sign  and  statistical  significance  for  the  in¬ 
dividual  models.  In  particular,  none  of  the  models  sim¬ 
ulates  an  increase  in  the  magnitude  of  the  trend  from 
1980-99  to  2000-49  and  from  2000-49  to  2050-99  that  is 
statistically  significant,  indicating  that  there  is  no  signifi¬ 
cant  increase  in  Arctic  PSC  occurrence  throughout  the 
twenty-first  century.  This  result  is  consistent  with  the 
small  changes  in  monthly  mean  temperatures  at  50  hPa 
that  are  simulated  by  the  models  in  January-February 
throughout  the  twenty-first  century  in  the  Arctic  (see 
Fig.  2). 

The  above  conclusions  are  based  on  daily  tempera¬ 
tures.  However,  while  the  temperature  for  NAT  and  STS 
formation  is  typically  about  195  K  in  the  stratosphere, 
the  exact  value  depends  on  the  partial  pressures  of  water 
vapor  and  nitric  acid  (Hanson  and  Mauersberger  1988). 


The  lower  panels  in  Fig.  4  show  the  actual  PSC  area  that 
is  calculated  by  the  model’s  PSC  schemes,  though  this 
was  only  available  from  five  models  (CCSR/NIES,  MRI, 
ULAQ,  UMSLIMCAT,  and  WACCM).  Both  in  the 
Antarctic  and  the  Arctic,  trends  in  the  actual  PSC  area 
simulated  by  these  models  are  similar  to  the  temperature- 
derived  PSC  trends,  except  for  the  UMSLIMCAT,  which 
shows  a  large  statistically  significant  positive  trend  be¬ 
tween  2000  and  2020,  but  a  negative  trend  in  Tnat.  The 
absolute  values  for  ?nat  (upper  panels  in  Fig.  4)  and 
the  actual  PSC  area  (lower  panels  in  Fig.  4)  differ  in  the 
models,  especially  in  the  UMSLIMCAT  and  the  WACCM 
where  the  actual  PSC  area  is  lower  than  the  temperature- 
derived  PSC  area.  In  addition,  the  magnitude  of  the  ac¬ 
tual  PSC  area  differs  among  the  models.  Table  4  shows 
that  there  is  a  good  correlation  between  the  trends  in 
Antarctic  PSC  areas  at  50  hPa  derived  from  daily  tem¬ 
peratures  below  195  K  (NAT)  and  the  actual  PSC  area. 
The  correlation  coefficients  are  calculated  for  the  five 
models  plus  the  multimodel  mean  for  the  periods  of 
1980-99,  2000-49,  and  the  entire  simulation,  and  from 
two  models  and  the  multimodel  mean  for  2050-99.  For 
the  five  models  that  provided  this  additional  PSC  diag¬ 
nostic  the  correlation  coefficient  between  the  Antarctic 
temperature-derived  and  the  actual  PSC  area  was  0.86, 
with  the  slope  of  the  least  squares  linear  fit  line  implying 
a  0.91%  increase  in  the  fraction  of  Tnat  for  each  percent 
increase  of  the  fraction  in  the  actual  PSC  area.  An  in¬ 
teresting  feature  is  that  the  correlation  between  7NAT  and 
the  actual  PSC  area  decreases  from  the  past  to  the  future, 
that  is,  the  correlation  coefficient  is  highest  for  1980-99, 
decreases  for  2000-49,  and  further  decreases  for  2050- 
99.  The  correlation  in  the  Arctic  between  trends  in  Tnat 
and  the  actual  PSC  area  is  smaller  than  the  correlation 
in  the  Antarctic  for  the  entire  simulation  but  is  higher 
for  the  shorter  periods  (1980-99,  2000-49,  2050-99;  see 
Table  4). 

c.  Zonal  mean  winds 

Changes  in  the  zonal  mean  eastward  wind  can  impact 
both  the  resolved  (Rossby)  wave  propagation  charac¬ 
teristics  in  the  models  and  the  filtering  of  the  para¬ 
meterized  subgrid-scale  gravity  waves  (see  section  4c). 
In  the  models  the  latitudinal  dependence  of  the  trends  in 
the  zonal  and  seasonal  mean  temperatures  (cf.  section 
3a  and  Figs,  lb  and  2)  alters  both  the  meridional  tem¬ 
perature  gradients  in  the  lower  stratosphere  and  the 
zonal  mean  zonal  winds,  thus  maintaining  thermal  wind 
balance.  Figure  6  shows  the  trends  in  the  eastward  winds 
at  70  hPa  for  the  two  solstice  seasons  for  those  simulations 
that  extended  to  at  least  2050.  Except  at  the  latitudes 
where  the  trends  change  sign,  all  of  the  trends  shown  are 
statistically  significantly  different  from  zero  at  the  95  % 
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Fig.  5.  (top)  Antarctic  (July-December)  and  (bottom)  Arctic  (November-April)  PSC  area 
at  50  hPa  (km2;  APSC)  from  the  models  (colors)  and  ERA-40  (black  dashed  line).  The  solid 
black  line  is  the  multimodel  mean  with  the  gray  shading  denoting  two  standard  errors.  The 
long-term  trend  in  the  simulations  and  the  ERA-40  data  is  calculated  from  the  maximum  values 
of  ApSC  for  5-yr  intervals,  which  are  displayed  as  circles.  Note  the  different  scale  between  the 
top  and  bottom  panels. 


confidence  level,  again  using  a  lag-1  autocorrelation  to 
account  for  the  nonindependence  of  the  residual  values 
about  the  trend  line. 

Although  there  is  quite  a  broad  spread  in  the  trends 
shown  in  Fig.  6,  most  of  the  simulations  exhibited  the 
same  basic  features  that  were  present  throughout  the 
lower  stratosphere  (e.g.,  similar  figures  were  obtained 


for  the  trends  at  100  and  50  hPa,  not  shown).  The  main 
exception  is  the  ULAQ  simulation,  especially  in  the 
Northern  Hemisphere.  Most  likely  this  is  a  consequence 
of  the  poor  representation  of  the  upper-tropospheric 
dynamics  in  that  model  resulting  from  its  quasigeostrophic 
formulation  and  coarse  vertical  resolution,  though  it 
should  also  be  noted  that  there  is  no  long-term  trend  in 
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TABLE  4.  Correlation  coefficient  and  slope  of  the  least  squares 
linear  fit  for  Tnat  and  the  actual  PSC  area. 


Trend  period 

Correlation/slope 
in  NH 

Correlation/slope 
in  SH 

Entire  simulation 

0.64/0.48 

0.86/0.91 

1980-99 

0.98/1.57 

0.90/1.03 

2000^19 

0.98/1.03 

0.85/0.59 

2050-99 

—0.67/— 0.36 

0.32/0.15 

the  SST  forcing  in  the  ULAQ  simulation  (see  Table  2). 
Because  of  these  deficiencies  the  results  from  the 
ULAQ  simulation  are  not  included  in  the  multimodel 
mean  shown  in  Fig.  6. 

In  the  subtropics  all  of  the  models  predict  an  eastward 
acceleration  in  the  zonal  mean  wind  in  both  hemi¬ 
spheres  and  seasons  corresponding  to  a  strengthening  of 
the  subtropical  jets  in  the  models.  For  the  Northern 
Hemisphere  winter  the  models  fall  into  two  categories. 
In  the  AMTRAC  and  MRI  simulations  the  eastward 
acceleration  extends  over  the  entire  extratropics  and 
peaks  in  the  midlatitudes.  In  the  CCSR/NIES,  CMAM, 
GEOS  CCM,  SOCOL,  and  WACCM  simulations  the 
eastward  acceleration  peaks  in  the  subtropics  and 


decreases  or  even  changes  sign  on  moving  toward  the 
pole.  The  reasons  why  the  models  fall  into  two  cate¬ 
gories  is  unclear,  but  most  likely  it  is  due  to  the  un¬ 
certainty  in  sign  of  the  high-latitude  response,  which  is 
strongly  influenced  by  changes  in  wave  driving  and  is 
not  robust  among  the  models.  In  contrast,  the  eastward 
acceleration  of  the  subtropical  jets  results  directly  from 
the  radiative-convective  warming  of  the  tropical  upper 
troposphere  and  therefore  consistently  has  the  same 
sign  for  all  of  the  models. 

In  the  high-latitude  Southern  Hemisphere  summer  all 
of  the  models  apart  from  the  SOCOL  model  indicate  a 
westward  acceleration  (Fig.  6a)  that  is  consistent  with  the 
warming  trends  in  the  lower  stratosphere  in  Fig.  2a,  re¬ 
sulting  from  the  recovery  in  the  Antarctic  ozone  [Fig.  3; 
see  also  Son  et  al.  (2008)].  Those  simulations  with  the 
strongest  warming  trends  in  Fig.  2  (e.g.,  WACCM  and 
GEOS  CCM)  were  also  the  simulations  with  the  largest 
westward  acceleration  apart  from  that  for  the  ULAQ 
model.  In  some  of  the  simulations  this  westward  trend 
in  the  seasonal  mean  zonal  winds  is,  in  part,  due  to  the 
earlier  transition  from  the  winter  eastward  jet  to  the 
summer  westward  jet.  This  is  consistent  with  other 
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Fig.  6.  Trends  in  the  zonal  mean  eastward  wind  at  70  hPa  as  a  function  of  latitude  for  (a)  December-February  and 
(b)  June-August.  The  trends  were  calculated  from  a  least  squares  linear  fit  to  the  modeled  time  series  of  the  zonal 
mean  eastward  wind  for  the  years  2000-99  for  the  AMTRAC,  CCSR/NIES,  CMAM,  GEOS  CCM,  and  MRI  sim¬ 
ulations  (solid  curves),  and  for  the  years  2000-50  for  the  SOCOL,  ULAQ,  and  WACCM3  simulations  (dashed 
curves).  The  extension  _e  to  a  model  short  name  indicates  that  the  ensemble  mean  of  more  than  one  simulation  was 
used  to  calculate  the  trend  for  that  model.  The  black  line  is  the  multimodel  mean  (excluding  ULAQ;  see  the  text), 
with  the  gray  shading  denoting  ±  two  standard  error. 
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results  that  have  shown  that  the  presence  of  an  Antarctic 
ozone  hole  can  delay  the  breakdown  of  the  Southern 
Hemisphere  polar  vortex  (e.g.,  Akiyoshi  et  al.  2009). 

4.  Brewer-Dobson  circulation 

a.  Tropical  upwelling 

Apart  from  the  short  (20  yr)  simulations  from  the  E39C 
and  MAECHAM4CHEM  models,  all  of  the  simulations 
show  a  clear  positive  trend  in  the  annual  mean  tropical 
upwelling  mass  flux  in  the  lower  stratosphere  (Fig.  7), 
consistent  with  previous  studies  (e.g.,  Butchart  and  Scaife 
2001;  Butchart  et  al.  2006;  Li  et  al.  2008;  McLandress  and 
Shepherd  2009).  As  in  Butchart  et  al.  (2006),  there  is  quite 
a  wide  spread  between  the  simulations  in  the  strength  of 
the  upwelling  (ranging  from  roughly  5  to  7  Mt  s~ 1  around 
the  year  20001  for  the  annual  mean;  see  also  Fig.  11), 
though  this  spread  is  much  smaller  (—0.2  Mt  s  ')  for 
simulations  from  the  same  model  (see  the  AMTRAC, 
CMAM,  E39C  model,  or  WACCM  results  in  Fig.  7, 
where  the  line  thickness  indicates  the  spread  between 
ensemble  members).  Indeed,  for  the  AMTRAC,  CMAM, 
E39C  model,  and  WACCM,  the  spread  between  ensem¬ 
ble  members  is  much  less  than  the  interannual  variability. 

For  the  majority  of  the  models  Fig.  7  shows  that  the 
amplitude  of  the  annual  cycle  is  larger  than  the  inter¬ 
annual  variability  in  either  the  seasonal  or  annual  mean 
upwelling.  In  general  the  strongest  upwelling  occurs  dur¬ 
ing  the  boreal  winter  (dark  gray  curves)  and  the  weakest 
during  the  boreal  summer  (light  gray  curves),  with  the 
amplitude  of  the  annual  cycle  typically  about  1  Mt  s  . 
The  SOCOL  and  ULAQ  results  were  an  exception,  with 
the  SOCOL  model  having  a  rather  weak  annual  cycle 
and  the  ULAQ  model  having  the  incorrect  phase.  The 
reasons  for  this  peculiar  behavior  are  unclear  because 
both  of  these  models  correctly  reproduced  the  annual 
cycle  in  the  100-hPa  temperatures  at  the  equator  (see 
Fig.  7  of  Eyring  et  al.  2006).  Figure  7  confirms  the  find¬ 
ings  of  Butchart  and  Scaife  (2001)  and  Butchart  et  al. 
(2006),  indicating  that  the  increases  in  tropical  upwelling 
occurs  throughout  the  year  (or  at  least  in  both  the  sol¬ 
stice  seasons  shown  in  the  figure). 

The  positive  trend  in  the  annual  mean  tropical  up¬ 
welling  is  a  robust  feature  of  the  multimodel  ensemble  at 
all  stratospheric  levels,  that  is,  the  predicted  strengthening 
of  the  Brewer-Dobson  circulation  is  not  a  shallow 


1  Note  that  the  AMTRAC  results  were  provided  at  77  hPa  and 
not  70  hPa,  and  therefore  the  upwelling  mass  flux  shown  for  the 
AMTRAC  in  Figs.  7  and  11  probably  slightly  overestimates  the 
mass  flux  at  70  hPa  because  the  upwelling  mass  flux  generally  de¬ 
creases  with  height  in  the  lower  stratosphere. 


feature,  but  extends  into  the  upper  stratosphere  (see  Fig. 
8a).  At  70  hPa  the  multimodel  mean  trend  is  just  under 
2%  decade and  slightly  less  than  that  estimated  by 
Butchart  et  al.  (2006),  though  the  difference  can  prob¬ 
ably  be  attributed  to  the  choice  of  year  with  respect  to 
which  the  percentages  were  calculated.  For  the  percent¬ 
age  trend  in  the  annual  mean  upwelling  (Fig.  8a)  all  of  the 
models  exhibit  the  same  features  with  a  minimum,  but 
still  positive,  trend  in  the  middle  stratosphere  and  an  in¬ 
creasing  trend  with  altitude  above  that.  Similar  behavior 
can  be  seen  in  the  trends  for  the  solstice  seasons  (Figs. 
8b,c)  though  there  is  a  greater  spread  between  the 
models,  and  one  model  (SOCOL)  predicts  a  negative 
trend  in  the  lower  to  middle  stratosphere  in  J une- August. 

b.  Age  of  air 

Mean  age  of  air  is  defined  as  the  mean  time  that  a 
stratospheric  air  mass  has  been  out  of  contact  with  the 
well-mixed  troposphere  and  provides  valuable  infor¬ 
mation  on  the  integrated  transport  within  the  strato¬ 
sphere.  For  the  REF2  simulations,  age-of-air  diagnostics 
were  available  in  the  central  database  for  the  following 
five  models:  AMTRAC,  CMAM,  GEOS  CCM,  ULAQ, 
and  WACCM.  When  these  models  were  used  to  simu¬ 
late  the  recent  past,  Eyring  et  al.  (2006)  found  that  at 
50  hPa  there  was  very  good  agreement  among  the 
models  and  with  the  age  inferred  from  observations  of 
trace  constituents.  However,  they  found  larger  differ¬ 
ences  in  the  upper  stratosphere,  both  between  the  models 
and  between  the  models  and  observations.  In  particular, 
the  AMTRAC,  CMAM,  GEOS  CCM,  and  WACCM  all 
had  a  lower  age  of  air  than  that  observed  and  the  ULAQ 
model  had  air  that  was  older  than  that  observed.  The 
mean  age  of  air  in  the  REF2  simulations  for  the  year 
2000  is  shown  in  Fig.  9  (see  figure  legend  for  details).  The 
overall  age  distributions  are  similar  in  all  models,  but 
there  are  quantitative  differences.  The  largest  deviations 
between  the  models  are  generally  in  regions  where  the  air 
is  really  old.  In  particular,  the  air  in  the  CMAM’s  and 
ULAQ’s  upper  stratosphere  and  mesosphere  and  in  the 
GEOS  CCM’s  extratropical  lower  stratosphere  is  more 
than  0.6  yr  older  than  the  multimodel  mean  age  (Fig.  9 
bottom  right)  of  4-5  yr  in  these  regions. 

Comparing  Figs.  7  and  9  suggests  that  to  a  large  extent 
the  intermodel  differences  in  the  mean  age  of  air  can 
be  attributed  to  the  differences  in  the  strength  of  the 
Brewer-Dobson  circulation  or  more  precisely  the  tropi¬ 
cal  upwelling  at  70  hPa.  The  two  models  with  the  stron¬ 
gest  upwelling  (the  AMTRAC  and  WACCM)  have  the 
youngest  air  (blue  shading  in  Figs.  9a,e).  For  the  five 
models  considered,  this  relationship  between  the  strength 
of  the  upwelling  and  mean  age  of  air  appears  to  be  true 
throughout  most  of  the  stratosphere  and  mesosphere,  but 
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Fig.  7.  Time  series  of  the  annual  (black  curves),  December-February  (dark  gray  curves),  and  June-August  (light  gray  curves)  mean 
tropical  upwelling  mass  flux  at  70  hPa  (77  hPa  for  the  AMTRAC  and  64  hPa  for  the  E39C  model)  for  each  model.  The  upward  mass  flux  is 
calculated  for  each  season  by  multiplying  the  seasonal  mean  residual  vertical  velocities  [W*;  see  Andrews  et  al.  1987,  their  Eq.  (3.5.1b)]  by 
a  scaling  proportional  to  the  cosine  of  latitude  then  integrating  over  all  latitudes  between  50°N  and  50°S,  where  W*  is  positive  (upward). 
The  density  at  70  hPa  was  assumed  to  be  70/1000  of  the  mean  surface  density  (1.2  kg  m~3).  The  annual  mean  mass  flux  is  the  average  of  the 
mass  flux  calculated  for  the  four  seasons,  and  therefore  takes  account  of  the  annual  cycle  in  which  the  upwelling  region  is  displaced  from 
the  equator  toward  the  summer  hemisphere.  The  thicker  curves  used  for  the  AMTRAC,  CMAM,  E39C,  and  WACCM  results  denote  the 
spread  between  the  ensemble  members.  (Note  the  mass  fluxes  were  only  available  for  two  of  the  three  WACCM  simulations.)  The  short 
20-yr  simulations  (E39C  and  MAECFIAM4CHEM)  and  also  those  of  the  SOCOL  and  ULAQ  models,  which  were  unable  to  reproduce 
the  correct  annual  cycle,  were  excluded  from  the  multimodel  mean  shown  in  the  bottom  middle  panel.  In  this  panel  the  shading  denotes  ± 
two  standard  error  for  the  annual  mean. 


5364 


JOURNAL  OF  CLIMATE 


Trend  in  upward  mass  flux  [%  per  decade]  Trend  in  upward  mass  flux  [%  per  decade] 


Volume  23 


-3-2-10  1  2  3 

Trend  in  upward  mass  flux  [%  per  decade] 


Fig.  8.  Vertical  profiles  of  the  trends  in  the  tropical  upwelling  for  models  with  simulations  longer  than  20  yr  for  the  (a)  annual  mean,  (b) 
December-February  mean,  and  (c)  June-August  mean.  At  each  pressure  level  the  upwelling  mass  flux  was  calculated  using  the  method 
described  in  Fig.  7.  The  trends  were  calculated  for  the  entire  time  series  shown  in  Fig.  7,  and  for  each  level  they  are  expressed  as 
a  percentage  of  the  upwelling  flux  obtained  from  the  least  squares  linear  fit  for  the  year  2000.  Levels  where  the  trends  are  statistically 
significant  at  the  95%  confidence  level  are  indicated  by  the  filled  circles.  The  black  curve  is  the  multimodel  mean  (excluding  the  SOCOL 
and  ULAQ  models;  see  the  caption  for  Fig.  7)  and  the  gray  shading  is  the  ±  two  standard  error. 


it  breaks  down  in  the  extratropical  lower  stratosphere. 
In  this  region  the  mean  age  of  air  is  determined  as  much 
by  the  horizontal  (isentropic)  mixing  as  by  the  large- 
scale  overturning  Brewer-Dobson  circulation.  There¬ 
fore,  the  older-than-average  air  in  the  GEOS  CCM’s 
extratropical  lower  stratosphere  is  probably  indicative  of 
weaker  mixing  in  that  model  compared  to  the  other  three 
models.  The  low  horizontal  resolution  of  the  CMAM 
and  the  WACCM,  and  particularly  the  ULAQ  model, 
likely  accounts  for  the  weaker-than-average  horizontal 
age  gradients  in  the  low-  to  midlatitude  lower  strato¬ 
sphere  in  these  three  models  [see  also  Fig.  10  of  Eyring 
et  al.  (2006)  for  the  CMAM  and  WACOM], 

In  all  five  models  the  annual  mean  of  age  of  air  de¬ 
creases  through  the  twenty-first  century  everywhere  in 
the  middle  atmosphere  (see  Fig.  10),  as  has  been  noted 
for  individual  models  (Austin  and  Li  2006;  Garcia  and 
Randel  2008;  Oman  et  al.  2009).  This  is  consistent  with 
an  acceleration  of  the  Brewer-Dobson  circulation  (see 
section  4a)  and,  for  the  AMTRAC,  in  agreement  with 
the  results  of  Austin  and  Li  (2006).  Also  in  all  of  the 
models  the  trends  are  statistically  different  from  zero  at 
the  95%  confidence  level. 

Above  the  lower  stratosphere  the  rates  of  decrease  in 
age  (negative  trend)  are  fairly  uniform  everywhere.  On 
average,  the  rate  of  decrease  is  much  less  (—20%)  in  the 
AMTRAC  than  in  the  other  models,  despite  the  fact 
that  there  is  little  difference  in  the  rate  of  increase  in 
tropical  upwelling  in  the  AMTRAC  compared  to  the 


CMAM,  GEOS  CCM,  or  WACCM  (see  Figs.  7  and  11). 
The  largest  rates  of  decrease  in  the  mean  age  are,  how¬ 
ever,  seen  in  the  lower  stratosphere  in  the  mid-  to  sub¬ 
tropical  latitudes,  especially  in  the  WACCM  and  in  the 
Northern  Hemisphere  in  the  GEOS  CCM.  This  suggests 
that  climate  change  is  producing  a  more  rapid  transport 
of  air  from  the  tropics  to  midlatitudes  in  the  lower 
stratosphere. 

While  the  models  all  show  a  decrease  of  age  of  air  with 
time,  for  simulations  of  both  the  recent  past  (not  shown) 
and  this  century  (see  also  Austin  and  Li  2006;  Garcia 
and  Randel  2008;  Oman  et  al.  2009;  Waugh  2009),  it  is 
important  to  note  that  a  recent  analysis  of  a  few  obser¬ 
vations  from  the  Northern  Hemisphere  midlatitude  lower 
stratosphere  (32°-51°N  and  24-35  km)  shows  a  weak  in¬ 
crease  in  estimates  of  age  of  air  (Engel  et  al.  2009).  How¬ 
ever,  the  large  uncertainties  in  this  observational  trend 
estimate  mean  that  the  observations  are  not  inconsistent 
with  the  model  trends  (Engel  et  al.  2009). 

c.  Wave  driving;  downward  control 

The  Brewer-Dobson  circulation  is  a  wave-driven 
circulation  (Holton  et  al.  1995).  For  steady  conditions, 
the  downward  control  principle  of  Haynes  et  al.  (1991) 
allows  the  residual  vertical  velocity  w*(<f>,  z)  at  latitude 
d>  and  log(pressure)-height  z  to  be  expressed  in  terms  of 
the  vertical  integral  of  the  mean  zonal  forces  T  acting 
above  z  as 
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where  m  =  acos<f>(u  +  «flcos</>)  is  the  angular  mo¬ 
mentum  per  unit  mass,  a  is  the  earth’s  radius,  fi  its  an¬ 
gular  velocity,  p0  a  basic-state  density,  and  T7  the  zonal 
mean  eastward  wind.  The  integration  is  up  a  line  of 
constant  zonal  mean  absolute  angular  momentum  but 
breaks  down  near  the  equator  where  the  lines  are  no 


longer  approximately  vertical  (Haynes  et  al.  1991).  None¬ 
theless,  the  total  upwelling  mass  flux  in  the  tropical  region 
can  be  derived,  by  continuity,  from  the  sum  of  the  ex- 
tratropical  downwelling  mass  fluxes  in  both  hemispheres 
that  are  calculated  using  this  downward  control  rela¬ 
tionship.  For  this  study  all  of  the  downward  control 
integrals  are  calculated  up  a  line  of  constant  latitude. 

In  the  middle  atmosphere  the  zonal  mean  forces  T  are 
due  to  Rossby  and  gravity  wave  breaking  and  other 
dissipative  eddy  processes.  In  the  models  the  resolved 
contribution  to  the  wave  breaking  is  proportional  to  the 
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Fig.  10.  Trend  in  the  annual  mean  age  of  air  for  the  AMTRAC,  CMAM,  GEOS  CCM,  ULAQ  model,  and 
WACCM  (yr  decade-1)  along  with  the  (bottom  right)  multimodel  mean  trend.  When  data  were  available  for  more 
than  one  integration  of  a  particular  model,  the  ensemble  mean  is  used  for  that  model. 


divergence  of  the  Eliassen-Palm  (EP)  flux  of  wave  ac¬ 
tivity  (e.g.,  chapter  3  in  Andrews  et  al.  1987).  The  con¬ 
tributions  to  T  from  subgrid-scale  gravity  waves  are 
parameterized  with  separate  schemes  for  the  waves  with 
orographic  and  nonorographic  sources  (see  section  2  and 
Table  2  for  details).  A  common  feature  of  all  of  the  pa- 
rameterizations,  however,  is  the  assumption  of  vertical 
wave  propagation.  Consequently,  for  the  parameterized 
subgrid-scale  waves,  the  downward  control  integral  sim¬ 
ply  depends  on  the  wave  momentum  flux  at  the  lower 
bound  of  the  integral  and  any  flux  that  is  allowed  to  es¬ 
cape  to  outer  space  (Shaw  and  Shepherd  2007)  and  does 
not  depend  on  details  of  where  the  momentum  is  de¬ 
posited  between  the  upper  and  lower  bounds.  In  gen¬ 
eral,  the  contributions  to  the  downward  control  integral 
fall  off  rapidly  with  height,  and  therefore  details  of  the 


momentum  flux  escaping  to  space  and  also  other  dissi¬ 
pative  processes,  such  as  Rayleigh’s  friction  or  diffusion 
acting  near  the  model’s  upper  boundary,  can  be  neglected 
when  calculating  w*  at  70  hPa  from  the  downward  con¬ 
trol  principle,  provided  that  the  model  has  a  sufficiently 
high  upper  boundary  and  momentum  conservation  is  not 
violated  within  the  stratosphere,  for  example,  by  failing 
to  deposit  orographic  GWD  momentum  (AMTRAC)  or 
employing  Rayleigh  friction  (ULAQ).  All  of  the  models 
used  in  this  study  had  an  upper  boundary  above  the 
stratopause  apart  from  the  E39C  model,  which,  for  sim¬ 
ulating  the  middle  atmosphere,  had  an  unusually  low 
upper  boundary  in  the  middle  stratosphere  (see  Table  2). 
Hence,  any  downward  control  calculations  shown  for 
the  E39C  model  at  70  hPa  have  to  be  treated  with  some 
caution.  Nonetheless,  Garny  et  al.  (2009)  showed  that 
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Fig.  11.  Comparison  of  the  (top)  annual,  (middle)  December-February,  and  (bottom)  June-August  mean  tropical  upwelling  mass 
fluxes  and  mass  flux  trends  at  70  hPa  (black  bars)  with  those  derived  from  downward  control  (gray  bars).  The  results  are  for  the  total 
upward  mass  flux  between  the  turn  around  latitudes  where  the  residual  vertical  velocity  changes  sign  from  upward  to  downward.  The  dark 
gray  shading  represents  the  contribution  from  the  EP  flux  divergence,  with  the  lighter  shadings  being  the  contributions  from  parame¬ 
terized  orographic  and  nonorographic  gravity  waves  (lightest  shade).  The  EP  flux  divergence  was  not  available  from  ULAQ;  GWD  was 
not  available  for,  or  not  used  (see  Table  1),  in  the  E39C,  SOCOL,  and  ULAQ  models;  in  addition,  nonorographic  GWD  was  unavailable 
for  the  AMTRAC.  The  GWD  results  shown  for  the  GEOS  CCM  are  the  combined  orographic  and  nonorographic  contributions.  The  solid 
horizontal  lines  are  the  multimodel  mean  with  the  dashed  lines  indicating  ±  two  standard  error.  When  results  were  available  for  more  than 
one  simulation  of  a  particular  model  the  ensemble  mean  is  used  in  the  figure.  The  diagonal  shading  indicated  negative  trends,  (left)  Mass 
fluxes  for  the  year  2000  based  on  a  linear  fit  to  the  time  series,  (right)  Trends  obtained  from  the  linear  fit. 
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analysis  of  the  mean  tropical  upwelling  mass  fluxes  in 
the  low  lid  model  E39C-A  at  76  hPa  derived  directly 
from  the  residual  velocity  and  calculated  using  down¬ 
ward  control  are  qualitatively  in  good  agreement. 

Figure  11  shows  that,  in  general,  there  was  reasonably 
good  agreement  between  the  actual  tropical  upwelling 
mass  fluxes  and  mass  flux  trends  and  those  derived  from 
the  downward  control  principle.  In  most  models  the 
orographic  gravity  waves  are  of  similar  importance  to 
the  resolved  waves  (EP  flux  divergence),  both  in  deter¬ 
mining  the  upwelling  and  its  trend.  For  the  five  models 
(AMTRAC,  CCSR/NIES,  CMAM,  GEOS  CCM,  and 
WACCM)  for  which  orographic  GWD  was  available 
above  70  hPa  (for  the  GEOS  CCM  the  total  GWD  was 
used)  this  accounted  for,  on  average,  30%  of  the  annual 
mean  tropical  upwelling  compared  to  67%  from  the  EP 
flux  divergence,  though  the  partitioning  varied  consid¬ 
erably  between  the  models  (from  14%  for  the  GWD  and 
88%  from  the  resolved  waves  in  the  GEOS  CCM  to  44% 
and  54%,  respectively,  in  the  CCSR/NIES  model).  In 
December-February  (June-August),  the  mean  contri¬ 
butions  from  the  orographic  gravity  waves  and  resolved 
waves  were  35%  (24%)  and  67%  (77%),  respectively. 
The  increased  contribution  from  orographic  GWD  dur¬ 
ing  the  boreal  winter/austral  summer  confirms  the  find¬ 
ings  of  Li  et  al.  (2008)  and  McLandress  and  Shepherd 
(2009)  and  establishes  that  this  aspect  is  robust  across 
a  range  of  models.  Unlike  the  orographic  GWD,  the 
nonorographic  GWD  had  very  little  impact  on  the  down¬ 
ward  control-derived  mass  fluxes  in  the  lower  stratosphere 
(see  light  gray  shading  in  Fig.  11)  because  most  of  the 
momentum  deposited  by  these  waves  in  the  extratropics 
occurs  well  above  the  70-hPa  level  (e.g.,  McLandress 
and  Shepherd  2009). 

As  first  noted  by  Li  et  al.  (2008)  for  the  AMTRAC 
simulations,  changes  in  orographic  GWD  can  be  a  sig¬ 
nificant  contributor  to  the  trends  in  the  upwelling,  espe¬ 
cially  during  the  boreal  winter.  The  contribution  from 
orographic  GWD  to  the  trend  in  upwelling  varied  con¬ 
siderably  over  the  five  models  with  GWD  diagnostics 
(AMTRAC,  CCSR/NIES,  CMAM,  GEOS  CCM,  and 
WACCM;  see  right-hand  panels  in  Fig.  11),  but  on  av¬ 
erage  it  accounted  for  59%  of  the  trend  in  the  annual 
mean  upwelling,  increasing  to  77%  during  December- 
February,  though  this  high  value  is  largely  due  to  just  two 
models  (the  AMTRAC  and  the  CCSR/NIES  model). 
Nonetheless,  even  the  smallest  23%  contribution  to  the 
increase  in  the  annual  mean  upwelling  from  GWD  in  the 
GEOS  CCM  is  large  when  it  is  considered  that  in  this 
model  only  a  modest  14%  of  the  driving  of  the  upwelling 
comes  from  GWD  (see  left-hand  panels  in  Fig.  11). 

Most  parameterizations  of  orographic  gravity  waves 
typically  deposit  momentum  in  the  westward  wind  shear 


on  the  upper  side  of  the  subtropical  jets,  especially  in  the 
Northern  Hemisphere.  This  region  includes  the  latitudes 
where  the  Brewer-Dobson  circulation  changes  from 
tropical  upwelling  to  extratropical  downwelling,  which 
are  the  latitudes  at  which  the  downward  control  in¬ 
tegrals  (see  equation  above)  were  calculated  to  derive 
the  strength  of  the  tropical  upwelling  shown  by  the  gray 
bars  in  Fig.  11.  Vertical  profiles  of  the  annual  mean 
orographic  GWD  at  these  “turn-around  latitudes”  (i.e., 
the  latitudes  where  Tv*  =  0)  are  shown  for  the  Northern 
Hemisphere  in  Fig.  12a  for  the  AMTRAC,  CCSR/NIES, 
CMAM,  GEOS  CCM,  and  WACCM  simulations  (see 
the  figure  legend  for  details).  At  the  Southern  Hemisphere 
turn-around  latitudes  the  vertical  profiles  are  similar,  but 
the  overall  orographic  GWD  is  much  weaker  and  is  not 
shown  here.  In  response  to  climate  change  in  the  simu¬ 
lations  there  was  an  eastward  acceleration  of  the  sub¬ 
tropical  jets  (see  section  3,  Fig.  6),  which  caused  the 
parameterized  gravity  waves  to  break  higher  in  the  at¬ 
mosphere  resulting  in  a  trend  for  more  drag  above  70  hPa 
and  less  below  as  shown  in  Fig.  12b.  In  turn,  the  increased 
orographic  GWD  above  70  hPa,  or  more  particularly 
the  increased  upward  gravity  wave  momentum  flux  at 
70  hPa,  contributed  significantly  to  the  strengthening  of 
the  tropical  upwelling  as  indicated  by  the  right-hand 
panels  in  Fig.  11.  Therefore,  in  the  case  of  the  vertically 
propagating  orographic  waves  it  is  the  secular  changes  to 
the  modeled  atmospheres  below  70  hPa  that  determined 
how  the  GWD  affected  the  response  of  the  tropical  up¬ 
welling  (at  70  hPa)  to  climate  change,  even  though  that 
response  was  the  result  of  “downward  control.” 

It  is  also  worth  noting  that  at  those  locations  where, 
and/or  times  when,  the  parameterized  orographic  waves 
are  already  breaking  (i.e.,  providing  drag)  below  70  hPa 
(cf.,  Fig.  12a  for  the  zonal  and  time  mean)  an  increase  in 
the  upward  momentum  flux  from  the  surface  cannot 
increase  the  flux  at  70  hPa  because  the  extra  flux  will  just 
contribute  to  the  momentum  deposited  at  the  level  at 
which  the  waves  first  break  or  saturate  (only  the  “un¬ 
saturated”  flux  is  allowed  to  propagate  upward  in  the 
parameterizations  schemes),  though  a  change  in  the  sur¬ 
face  flow  could  change  the  spatial  and  temporal  extent 
of  the  sources.  For  the  CMAM  there  was  only  a  small 
change  in  the  momentum  flux  from  the  surface  (not 
shown)  and  the  increased  drag  above  70  hPa  in  that 
model  came  mainly  from  reduced  filtering  (breaking) 
below  70  hPa  (McLandress  and  Shepherd  2009).  In  the 
WACCM  the  flux  from  the  surface  was  constant  and  for 
the  other  models  it  was  not  diagnosed. 

In  all  of  the  simulations  climate  change  also  increased 
the  EP  flux  entering  the  stratosphere  at  70  hPa  in  the 
subtropics  (Fig.  13),  consistent  with  an  increase  in  the  EP 
flux  divergence  above  that  level  accelerating  the  tropical 
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Fig.  12.  (a)  For  the  simulations  of  the  AMTRAC,  CCSR/NIES  model,  CMAM,  GEOS  CCM,  and  WACCM,  the 
vertical  profiles  of  the  annual  mean  orographic  gravity  wave  drag  [m  s-1  (day)  ']  at  the  latitudes  north  of  the 
equator,  where  W*  =  0  at  70  hPa.  The  annual  mean  is  the  average  of  the  four  seasonal  means  and  therefore  takes 
account  of  the  seasonal  cycle  in  the  position  where  W*  =  0.  The  results  shown  are  for  the  year  2000  based  on  a  linear 
fit  to  the  time  series  of  annual  mean  vertical  profiles,  (b)  Change  in  drag  (%  decade-1)  based  on  the  linear  fit  to  the 
time  series.  When  results  were  available  from  more  than  one  simulation  of  a  particular  model,  the  ensemble  mean  is 
shown.  In  the  AMTRAC  orographic  gravity  wave  drag  was  set  to  zero  above  42  hPa.  The  curves  for  the  AMTRAC 
and  the  CMAM  do  not  extend  to  the  lowest  levels  because  data  were  not  provided  for  those  levels.  Nonetheless,  it 
should  be  emphasized  that  the  orographic  GWD  does  operate  throughout  the  depth  of  the  troposphere  in  these 
models. 


upwelling  (e.g.,  dark  gray  bars  in  the  right-hand  panels  of 
Fig.  11).  However,  because  the  resolved  waves  propagate 
both  horizontally  and  vertically  and  can  be  refracted, 
it  was  not  possible  to  quantify  the  contribution  to  the 
downward  control  integral  in  terms  of  just  this  upward 
wave  momentum  flux  (i.e.,  vertical  component  of  the  EP 
flux)  at  70  hPa,  as  was  done  for  the  parameterized  waves. 
Similarly,  the  possibility  of  refraction  and  the  more  com¬ 
plex  nature  of  the  sources  of  the  resolved  waves  (Nathan 
and  Cordero  2007),  compared  to  those  used  for  the  pa¬ 
rameterized  waves,  make  it  difficult  to  identify  the  causes 
of  the  trend  in  upward  wave  flux  at  70  hPa.  Nonetheless, 
a  rather  detailed  analysis  of  the  CMAM  by  McLandress 
and  Shepherd  (2009)  showed  that  the  large-scale  waves 
(zonal  harmonics  1-3)  accounted  for  about  60%  of  the 
driving  of  the  annual  mean  upwelling  trend  by  resolved 
waves,  with  the  remainder  coming  from  synoptic  waves 
(zonal  harmonics  >3).  This  contribution  from  the  syn¬ 
optic  waves  increased  to  around  50%  during  December 
through  February.  Using  additional  runs  of  the  GEOS 
CCM  Oman  et  al.  (2009)  identified  increases  in  tropical 
SSTs  as  a  significant  contributor  to  an  increase  in  the 


resolved  wave  fluxes,  which  is  consistent  with  the  result 
of  Fomichev  et  al.  (2007)  that  the  associated  increase 
in  SSTs  adds  an  additional  1-K  cooling,  from  upwell¬ 
ing,  to  the  cooling  of  the  tropical  lower  stratosphere 
from  a  doubling  of  C02.  Neither  of  these  studies  at¬ 
tempted  to  separate  out  the  direct  contribution  resulting 
from  a  change  in  the  source  of  waves  from  the  indirect 
contribution  resulting  from  a  change  of  propagation  due 
to  the  tropospheric  climate  change.  Nonetheless,  the 
results  of  Oman  et  al.  (2009)  are  consistent  with  other 
studies  (e.g.,  Hardiman  et  al.  2007)  that  have  shown 
a  correlation  between  variability  in  tropical  SSTs  and 
tropical  upwelling  at  different  levels  though,  critically, 
Hardiman  et  al.  (2007)  noted  that  the  sign  of  the  cor¬ 
relation  changes  on  moving  from  the  troposphere  to  the 
stratosphere.  Deckert  and  Dameris  (2008)  have  also 
reported  that  higher  tropical  SSTs  can  lead  to  a  localized 
seasonal  strengthening  of  upwelling  in  the  tropics,  but  it 
was  unclear  from  the  results  they  show  if  this  would 
produce  a  net  increase  in  total  upwelling  between  the 
turn-around  latitudes  at  70  hPa  as  noted  in  sections  4a 
and  4c. 
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Fig.  13.  Time  series  of  the  annual  mean  upward  component  (F-)  of  the  EP  flux  at  70  hPa 
averaged  between  20°-40°N  and  20°-40°S  for  the  simulations  from  the  AMTRAC,  CCSR/NIES 
model,  CMAM,  GEOS  CCM,  MRI  model,  and  WACOM.  When  result  were  available  from 
more  than  one  simulation  of  a  particular  model  the  ensemble  mean  is  shown. 


5.  Discussion  and  conclusions 

In  this  study  the  response  of  the  stratospheric  climate 
and  circulation  to  increasing  amounts  of  GHGs  and  ozone 
recovery  in  the  twenty-first  century  has  been  assessed 
using  results  from  an  ensemble  of  near-identical  transient 
simulations  of  11  chemistry-climate  models  (CCMs).  In 
all  of  the  simulations  the  abundance  of  GHGs  increased 
according  to  the  IPCC  SRES  A1B  scenario  (Nakicenovic 
and  Swart  2000),  while  the  ozone  recovery  mainly  re¬ 
sulted  from  a  reduction  in  the  amounts  of  ozone-depleting 
substances  based  on  the  WMO  (2003)  Ab  scenario  for 
halogens,  though,  in  each  simulation,  the  recovery  of 
ozone  was  modulated  by  concomitant  changes  in  the 
modeled  climate  and  circulation  (see  Eyring  et  al.  2007). 
Consequently,  a  significant  difference  in  the  radiative 
forcing  of  the  climate  in  the  simulations  arose  from  the 
differences  in  the  model  projections  of  ozone  recovery 
that  were  reported  in  Eyring  et  al.  (2007).  Another 
important  model-dependent  difference  in  the  radiative 
forcing  resulted  from  the  sensitivity  of  stratospheric 
water  vapor  to  the  modeled  tropical  tropopause  temper¬ 
atures.  In  addition,  there  were  differences  in  the  forcing 
from  the  prescribed  SSTs  because  of  different  climate 
sensitivities  in  the  coupled  models  used  to  generate 
those  SSTs,  though,  again,  for  each  simulation  these 
evolved  consistently  with  the  SRES  A1B  scenario  for 
the  GHGs  (i.e.,  they  were  taken  from  coupled  ocean- 
atmosphere  model  projections  based  on  the  same  GHG 
scenario).  All  of  the  simulations  included  these  anthro¬ 
pogenic  forcings,  while  for  most  of  them  natural  forcings 
were  generally  excluded.  However  simulations  from 
three  of  the  models  included  a  future  solar  forcing  and 
two  models  included  an  artificial  QBO  forcing  in  the 
tropics  (see  columns  4  and  5  in  Table  2).  For  these  models 
and  simulations  it  was  assumed  that  any  nonlinear  recti¬ 
fication  of  this  artificially  forced  natural  variability  did  not 


contribute  directly  to  the  long-term  secular  changes,  and 
hence  the  changes  or  trends  reported  here  represent 
the  stratospheric  response  to  the  changes  in  anthropo¬ 
genic  forcings  (i.e.,  changes  in  the  amounts  of  GHGs  and 
ODSs). 

A  rather  important  component  of  the  “fingerprint”  of 
global  climate  change  is  a  cooling  of  the  stratosphere, 
which  increases  with  height.  While  ozone  depletion  has 
been  a  significant  contributor  to  global  mean  stratospheric 
cooling  in  the  past,  future  stratospheric  global  temperature 
changes  will  be  driven  overwhelmingly  by  climate  change 
(Shepherd  and  Jonsson  2008).  All  of  the  simulations 
assessed  in  this  study  produced  this  fingerprint  for  the 
global  temperature  trend  with  very  little  spread  between 
the  models  or  simulations.  On  average  the  annual  mean 
cooling  at  10  hPa  was  0.59  K  decade-1  with  the  standard 
error  of  the  ensemble  mean  of  0.037  K  decade-1.  The 
simulations  also  suggest  additional  spatial  and  temporal 
structure  to  the  trend,  which,  in  particular,  manifests  as 
a  warming  of  the  polar  lower  stratosphere  in  the  austral 
summer  with  enhanced  cooling  above  and  a  near-zero 
temperature  trend  in  the  polar  lower  stratosphere  in  the 
boreal  winter.  This  spatial  and  temporal  structure  re¬ 
sulted,  however,  from  a  combination  of  both  climate 
change  and  ozone  recovery. 

The  warming  over  Antarctica  has  been  noted  before  for 
individual  models  (e.g.,  Perlwitz  et  al.  2008;  McLandress 
and  Shepherd  2009)  and  for  IPCC  AR4  runs  with  pre¬ 
scribed  ozone  changes  (Son  et  al.  2008).  though  in¬ 
terestingly  the  signal  is  not  as  strong  in  the  AR4  runs  as 
in  the  chemistry-climate  model  simulations  with  their 
higher  upper  boundaries  (see  also  Shaw  et  al.  2009).  The 
results  presented  in  this  paper  confirm  this  to  be  a  robust 
feature  of  stratospheric  change  over  the  next  50  years 
or  so.  Moreover,  the  correlation  (0.83  for  the  period  of 
2000-50)  between  the  rate  of  recovery  of  the  Antarctic 
ozone  hole  (e.g.,  Eyring  et  al.  2007)  and  the  strength  of 
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the  polar  warming  in  the  simulations  implies  that  the 
latter  is  more  or  less  a  direct  radiative  response  to  ozone 
recovery,  although  it  is  also  important  to  look  at  the 
time  signature  of  the  changes  to  confirm  the  attribution 
(Hitchcock  et  al.  2009).  At  same  time  it  also  indicates 
a  sensitivity  of  the  temperature  response  to  details  and 
the  accuracy  of  the  interactive  ozone  chemistry  and  as¬ 
sociated  transport  of  constituents  in  the  models. 

Unlike  for  the  Antarctic,  there  has  been  little  con¬ 
sensus  on  how  the  Arctic  stratosphere  will  respond  to 
climate  change  and  ozone  recovery  (e.g.,  Austin  et  al. 
2003;  WMO  2007,  chapter  5).  However,  the  results  from 
the  longer  simulations  assessed  here  suggest  that  the 
extra  radiative  cooling  from  growing  amounts  of  GHGs 
will  be  approximately  balanced  by  a  concomitant  in¬ 
crease  in  the  adiabatic  warming  through  increased  polar 
downwelling.  The  net  effect  is  a  near-zero  temperature 
trend  in  the  Arctic  winter  lower  stratosphere.  Because 
this  signal  of  climate  change  is  the  small  residual  of  much 
larger,  but  similarly  sized,  perturbations  to  the  adiabatic 
warming  and  diabatic  cooling,  it  is  not  statistically  sig¬ 
nificant  in  a  single  simulation  with  respect  to  the  large 
natural  interannual  variability  in  this  region.  Nonetheless, 
the  results  from  the  multimodel  ensemble  confirm  that  the 
signal  is  probably  a  robust  feature  of  climate  change. 

An  important  consequence  of  the  long-term  strato¬ 
spheric  temperature  changes  in  the  polar  latitudes  is  the 
effect  on  the  occurrence  of  PSCs,  and  hence  ozone  loss 
in  these  regions.  In  this  study  trends  in  the  occurrences 
of  PSCs  inferred  from  the  modeled  daily  temperature 
fields  generally  agreed  well  with  the  trends  obtained 
directly  from  the  parameterizations  of  PSC  processes 
within  the  models.  A  fairly  robust  feature  in  the  Antarctic 
was  a  steep  increase  in  the  accumulated  area  of  PSCs  in 
response  to  the  development  of  the  ozone  hole,  with 
most  of  this  increase  consequently  coming  in  late  spring 
(Hitchcock  et  al.  2009),  followed  by  a  gentle  increase 
in  the  future  as  cooling  from  climate  change  dominates 
the  warming  from  ozone  recovery.  Because  the  cli¬ 
mate  change-induced  changes  in  PSC  occurrence  oc¬ 
cur  throughout  the  winter  (Hitchcock  et  al.  2009),  this 
means  that  even  as  the  ODSs  begin  to  disappear  from 
the  stratosphere  the  areal  extent  of  the  ozone  hole  could 
still  expand  further  before  it  disappears.  In  the  Arctic 
there  was  much  greater  spread  between  the  models,  but 
consistently  with  the  near-zero  seasonal  mean  tempera¬ 
ture  trends  in  this  region  (see  above),  the  models  suggest 
that  on  average  there  would  be  no  significant  change  in 
occurrences  of  PSCs  in  the  near  future  and  throughout 
the  twenty-first  century. 

Another  important  consequence  of  the  temperature 
trends  was  a  change  in  the  temperature  gradients  and 
thus,  through  thermal  wind  balance,  a  change  in  the  zonal 


mean  eastward  wind.  Results  from  the  multimodel  com¬ 
parison  indicate  two  robust  features.  First,  there  was 
an  eastward  acceleration  of  the  subtropical  jet  in  both 
hemispheres,  which  had  important  ramifications  for  both 
the  parameterized  gravity  waves  and  also  the  resolved 
(Rossby)  waves  propagating  into  the  stratosphere  (see 
below).  Second,  there  was  a  westward  acceleration  of 
the  lower-stratospheric  winds  over  Antarctica  during 
December-February,  which  corresponds  to  a  decrease 
in  the  amplitude  of  the  southern  annular  mode.  Like  the 
temperatures  in  this  region  this  response  was  sensitive  to 
the  strength  of  the  ozone  recovery  in  the  models  and  was 
not  present  in  those  twenty-first  century  IPCC  AR4  sim¬ 
ulations  that  did  not  include  the  ozone  recovery  or  had 
a  poor  representation  of  the  stratosphere  (Perlwitz  et  al. 
2008;  Son  et  al.  2008).  This  is  particularly  significant 
because  there  is  evidence  (Perlwitz  et  al.  2008;  Son  et  al. 
2008)  that  indicates  the  signal  in  the  stratospheric  winds 
can  extend  to  the  surface  where  there  are  possible  inter¬ 
esting  implications  for  the  carbon  uptake  by  the  Southern 
Ocean  (e.g.,  Matear  and  Lenton  2008). 

In  all  of  the  models  the  annual  mean  tropical  up- 
welling  at  70  hPa  was  predicted  to  increase  by,  on  av¬ 
erage,  —2%  decade-1,  confirming  the  earlier  results  of 
Butchart  et  al.  (2006).  The  results  of  this  study  also  es¬ 
tablished  that  the  predicted  increase  occurred  throughout 
the  depth  of  the  stratosphere  and  caused  a  fairly  uniform 
reduction  in  the  age  of  air  of  roughly  0.05  yr  decade-1 
in  those  models  for  which  this  diagnostic  was  available. 
Moreover,  the  differences  between  the  mean  age  in 
these  models  correlated  well  with  respective  differ¬ 
ences  in  the  tropical  upwclling  (i.e.,  models  with  a 
faster-than-average  Brewer-Dobson  circulation  had 
younger-than-average  air  and  vice  versa),  apart  from  in 
the  extratropical  lower  stratosphere  where  horizontal 
mixing  is  as  important  as  the  large-scale  overturning 
circulation.  This  suggests  that  the  spread  in  the  simu¬ 
lated  age  can  be  attributed  mostly  to  a  corresponding 
spread  in  the  strength  of  the  Brewer-Dobson  circulation 
and  not  directly  to  the  differences  in  the  tracer  advection 
schemes.  On  the  other  hand,  the  differences  in  the  ad¬ 
vection  schemes  could  contribute  to  the  differences  in 
horizontal  mixing. 

In  those  models  for  which  the  information  was  avail¬ 
able  30%  of  the  driving  and  59%  of  the  trend  in  the  an¬ 
nual  mean  tropical  upwelling  in  the  lower  stratosphere 
was,  on  average,  due  to  the  parameterized  orographic 
gravity  waves.  The  importance  of  the  orographic  as  op¬ 
posed  to  nonorographic  gravity  waves  was  a  common 
feature  of  all  the  models  that  provided  these  diagnostics. 
While  it  should  be  noted  that  the  momentum  fluxes  in 
the  lower  stratosphere  resulting  from  orographic  gravity 
waves  are  not  well  constrained  by  observations,  the  basic 
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mechanism  of  orographic  gravity  wave  breaking  heights 
shifting  upward  as  the  zonal  winds  strengthen  is  based  on 
fundamental  fluid  dynamical  principles.  With  the  for¬ 
mulations  used,  the  trend  in  upwelling  resulting  from  the 
orographic  waves  is  a  direct  consequence  of  the  eastward 
acceleration  of  the  subtropical  jets  in  response  to  climate 
change  (see  above)  allowing  more  waves  to  propagate 
and  break  above  70  hPa. 

The  overall  conclusion  of  this  analysis  is  that  there  are 
several  robust  features  that  all  (or  nearly  all)  models 
predict,  that  is,  consistent  trends  in  temperature,  PSC 
occurrence,  zonal  winds,  tropical  upwelling,  age  of  air, 
and  the  processes  driving  the  Brewer-Dobson  circula¬ 
tion.  Further  improvement  of  the  models,  including  their 
representation  of  tropospheric  climate  change,  should 
lead  to  a  refinement  of  these  predictions.  However,  the 
trends  in  some  quantities  (e.g.,  lower-stratosphere  south 
polar  temperatures)  are  sensitive  to  the  projected  rate  of 
ozone  recovery  in  the  models,  which  depends  on  the 
extent  of  their  simulated  ozone  depletion.  Hence,  im¬ 
proving  the  predictions  of  these  quantities  will  depend 
on  improving  the  simulation  of  ozone  depletion  and 
recovery  and  will  require  the  additional  consideration 
of  the  chemistry  and  transport  characteristics  of  the 
models,  which  is  quite  beyond  the  scope  of  this  study. 
Furthermore,  to  properly  assess  the  impact  of  the  changes 
in  the  stratosphere  on  the  surface  climate  will  require 
simulations  with  a  fully  interactive  ocean. 
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